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Using photon conversions to e + e~ pairs, the energy spectrum of inclusive photons from ip(2S) 
radiative decays is measured by BESII at the Beijing Electron-Positron Collider. The Xcj(lf) 
states (J=0,l,2) are clearly observed with energy resolution between 2.3 to 3.8 MeV, and their 
masses and the spin-averaged Xcj mass are determined to be M Xc0 — 3414.21 ± 0.39 ± 0.27, M Xal = 
3510.30 ± 0.14 ± 0.16, M Xc2 = 3555.70 ± 0.59 ± 0.39 and M( 3 P cog ) = 3524.85 ± 0.32 ± 0.30 MeV/c 2 . 
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PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx, 13.40Hg 

I. INTRODUCTION 

Precise measurements of the spectrum and the decay properties of charmonia are essential to test Potential QCD 
models and QCD based approaches Q. There is renewed interest since the discovery of the X(3872) Q and the 
observations of the expected r]c{2S) and he Q~Pi) states Q, and there has been recent progress, both theoretically 
and experimentally 4]. There are more accurate determinations of the charmonium mass spectrum and radiative 
transition rates using both a relativistic quark model with relativistic corrections of order v 2 /c 2 [jj and a potential 
model with a semirelativistic approach |6J]. The ip(2S) mass and width have been redetermined with an updated 
radiative correction 0, and newly measured with better precision @. In addition to previous measurements of 
Xcj states 0, two XcO measurements by E835 0| an d new Xcj (J=0,l,2) measurements by CLEO have been 
recently published. Improved precision on \cj masses is important for the determination of the singlet-triplet splitting, 
Af( 1 F 1 ) - M( 3 P cog ), which is predicted by lattice QCD and nonrelativistic QCD Here M( 3 P cog ) is the spin- 

averaged 3 Pj mass for the Xcj states (J = 0, 1, 2). 

In this paper, results on the Xcj masses (J=0,l,2) and widths (J=0,1) from a measurement of the energy spectrum of 
inclusive photons in ^(2S) radiative decays, using photon conversions to improve the energy resolution, are presented. 
The measurement uses 14 x 10 6 ip(2S) events collected with the upgraded Beijing Spectrometer (BESII) at the BEPC 
Collider. 



II. BES DETECTOR AND MONTE CARLO SIMULATION 

The BESII detector is described elsewhere ^3j. A 12-layer vertex chamber (VC) surrounding the beam pipe provides 
hit information in trigger criteria for charged tracks. Charged particle momenta are determined with a resolution 
of a p /p = 1.78% x /T + p2 ( p in QeV/c) in a 40-layer cylindrical drift chamber (MDC). Particle identification is 
accomplished by measurements of ionization (dE/dx) in the MDC and time-of-flight (TOF) in a barrel- like array of 
48 scintillation counters. The dE/dx resolution is o dE j dx = 8%; the TOF resolution is cttof=200 ps for hadrons. A 

12-radiation-length barrel shower counter (BSC) measures energies of photons with a resolution of cte/ 1 E — 21%/ \fE 
(E in GeV). A solenoidal coil supplies a 0.4 Tesla magnetic field over the tracking volume. 

A Geant3 based Monte Carlo (MC) SIMBES |l4|. which simulates the detector response including interactions 
of secondary particles in the detector material, is used to determine the energy resolution and detection efficiency 
of photons reconstructed from their converted e + e~ pairs, as well as to optimize selection criteria and estimate 
backgrounds. Under the assumption of a pure El transition for the ip(2S) — > ^XcJ decays, the polar angle (6) 
distributions of the photons are given by 1 + fccos 2 9 with k = 1,—^, j§ for J = 0, 1, 2, respectively [l5| |. 

Good energy resolution for low energy photons is essential for precise measurements of Xcj masses and widths 
from fitting the photon spectrum of tp(2S) radiative decays. Momentum resolution of about 1.6 to 4.1 MeV/c can 
be obtained for low momentum electrons from 60 to 250 MeV/c. Photons from ip(2S) — > "fXcj decays have energies 
of about 261, 171, and 128 MeV for the Xcj final states (J=0,l,2), and the electrons produced in photon conversions 
occur in this low momentum region. 

III. PHOTON RECONSTRUCTION AND SELECTION 

We choose two oppositely charged tracks with each track having a good helix fit and a polar angle with | cos 9\ < 0.8. 
The intersection of the electron and positron trajectories in the xy-plane (the beam line is the z axis) is determined, 
and this point is taken as the photon conversion point (CP). The photon conversion length R xy is defined as the 
distance from the beam line to the CP in the xy-plane. Fig. ^ shows the R xy distribution for photon conversions 
to e + e _ pairs in the BESII detector for hadronic events in the 58 xlO 6 J/ip event sample. The two broad peaks 
in Fig. ^ correspond to the beampipe region, where the beampipe, the VC, and inner wall of the MDC are located. 
Combinatorial background from charged hadron tracks is also seen in the R xy < 2 cm region. Equivalent materials 
in the beampipe wall, VC, VC outer wall, and the MDC inner wall are 0.536, 0.657,0.375, and 1.107 in units of 
O.OlAo 0], respectively, where Xq is a radiation length. The electron and positron directions are calculated at the 
photon conversion point, and their momenta are corrected to that point. 

Good photons are selected. The photon conversion length must lie within the beampipe region, 2 < R xy < 22 cm, 
and the invariant mass of an e + e~ pair is required to satisfy M e + e - < 20 MeV/c 2 . Combinatorial background from 
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FIG. 1: The R xy distribution for gamma conversions to e + e pairs in the BESII detector from hadron events in the 58 xlO 6 
J /i\) event sample. 

charged hadron tracks is further removed by requiring cosOdefi > 0.9 , where Odefi is the deflection angle between the 
photon momentum and photon track (a vector from the beam to the CP). To suppress background from beam-gas 
and beam-pipe interactions, the total energy in the event must satisfy E to t > Eb eam /2 and momentum asymmetry 
must satisfy dp asym < 0.9. Here dp asym is defined as a ratio of the vector sum to the scalar sum of the momenta 
of all charged and neutral tracks in the event. The observed photon energy spectrum from the ^(2S I ) data after the 
selection of good photons is shown in Fig. [3 The spectrum shows the XcJ states plus a large background. The sharp 
drop at low energy is mainly caused by low photon detection efficiency. 

.—. 1200, — , — , — , — , — , — , — , — , — , — , — , — , — , — , — | 
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FIG. 2: Inclusive photon spectrum from photon conversions from 14 x 10 6 ip(2S) events. A fit (described in the text) is made 
to ip(2S) — > 7Xcj decays (J = 0, 1, 2) plus threshold background. Points with error bars are data. The solid line is the fit; the 
dashed line is the fitted background. Background subtracted results are shown in the lower plot. 



IV. dE/dx CORRECTION AND PHOTON ENERGY SCALE 

Energy loss dE/dx by ionization for electrons traversing a small thickness of material with energy above a few tens 
of MeV can be described by the Bethc-Bloch equation :9j. The dE/dx correction for charged particles, produced near 
the beamline and traversing the whole beampipe region, should take into account the full thickness of material in the 
region. However, e + e~ pairs from photon conversions are mostly produced in the region where the VC outer wall 
and the MDC inner wall are located. Thus the effective thickness of material between the location, where a pair is 
produced, and first layer of the MDC wires must be estimated for each electron pair. The procedure to make dE/dx 
corrections for electrons has two steps : (1) A preliminary dE/dx correction using half the full thickness of all the 
materials in the beampipe region is made. Good photons are reconstructed, and their conversion lengths R xy are 
calculated. (2) The final dE/dx corrections are estimated based on these R xy . 

The energy scale of photons reconstructed from e + e~ pairs is studied using simulated MC events and data. 63 x 10 6 
7T° signal events are generated using MC technique, with same momentum and polar angular distributions as that 
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found from w° data sample. A sample of 7r° mesons decaying to two photons with both photons converting to e + e~ 
pairs is selected from 58 x 10 6 J/ip events. To suppress hadron contamination, electron identification is required 
and good photons are selected. Background is further suppressed with additional requirements on the photon energy, 
E~f < lGeV, and the opening angle between the two photons, 0.75 < | cos# 77 | < 0.97. The invariant mass distribution 
of two photons for both MC and data, after the specific dE/dx correction for electrons described above, is fitted with 
the improved Crystal Ball (ICB) function (defined in section^) plus a first order polynomial background. The results 
are shown in Fig. El The resulting tt° masses (134.47 ± 0.42) MeV/c 2 in data and (134.86 ± 0.20) MeV/c 2 in MC are 
consistent with the PDG value of 134.98 MeV/c 2 The corresponding mass resolutions (5.70 ± 0.58) MeV/c 2 in 
data and (5.55 ± 0.21) MeV/c 2 in MC agree within errors. The x 2 /D-F. (degree of freedom) for the fits are 126/103 
in data and 117/140 in MC. 

i 1 1 1 i 1 1 1 i 1 

i\ -_ 
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FIG. 3: Invariant mass distribution of photon pair from J/ip data (top) and Monte Carlo (bottom) events. The solid line is 
the fitted curve for signal plus background. The dashed line is the fitted curve for background. 
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V. PHYSICS PROBLEM AND DETECTOR RESOLUTION 



The energy of a photon from ip(2S) — > "fXcj decay is given by 

= (Af 2 (2S) -M 2 c; )/2A^ (2S) , (1) 

where M^ps) and M XcJ are the masses of the ip(2S) and Xcj, respectively. The i/j(2S) and \ c j widths must be taken 
into account; M^s) ancl M XcJ are described by Breit-Wigner functions (2D problem). By taking x = M^^sjt the- 
probability density function (pdf) for the photon energy E 7 can be written as |l7| 

f P df(E,) = J BW{x)BW{M XcJ )^—dx, (2) 

where M XcJ depends on E 1 by Eq. 0). 

As a result of traversing material in the beam pipe region, the electron energy is smeared due to energy loss 
by ionization, and a long tail on the low side of the energy distribution is induced by bremsstrahlung radiation. 
Multiple scattering of electrons, especially at large angles, gives tails on both sides of the photon energy distribution 
of photon conversions. The photon energy resolution from photon conversions can be nicely modeled by our Geant3 
MC simulation, and well fitted by the ICB function. The original Crystal Ball (CB) function has a Gaussian in its 
central and upper energy region but long tail at lower energy region |18|. The improved CB function is defined as 
same as the CB function but has an additional tail at its upper side. The photon energy distributions from large MC 
samples of ip(2S) — > jXcJ decays (J = 0, 1, 2), with zero widths for both the ip{2S) and %cj states, are fitted to ICB 
functions and shown in Fig. The x 2 /D.F. from the fits are 103.8/93, 37.7/53 and 47.6/43 corresponding to XcO, 
Xcij Xc2 states. Five parameters in the ICB function, the photon energy resolution and four empirical parameters 
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to describe the tails on the lower and upper sides are determined from the hts and used as input parameters in the 
detector resolution function for each decay mode. Photon energy resolutions for the ijj(2S) — > 7Xej decays (J=0,l,2) 
are found to be 3.78 ± 0.04, 2.58 ± 0.05, and 2.26 ± 0.11 MeV, respectively. 
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FIG. 4: Energy distributions of signal photons from t/>(2S) — > 7Xcj decays with zero widths of both ip(2S) and \ c j for Xco(left), 
Xci(mid) and Xc2 (right) final states fitted to the ICB function. The points are MC data. The solid line are the fitted curve. 

The energy dependencies of the photon detection efficiency and resolution are determined using MC simulation in 
the energy range between 100 and 400 (300, 220) MeV for ip(2S) — > 7Xc0 (lXci,lXc2) decays. The efficiency includes 
the effects of detector geometry, MDC tracking, photon reconstruction, and the spin-dependent cos 8 distribution. 



VI. FITTING THREE PHOTON SPECTRUM 



Background contamination from Xco — * 7 'J/ip decay is negligible due to its small branching fraction. To exclude 
photons from the Xci,c2 — > jJ/ip decays, the photon energy range in the fits is chosen to be within 90 MeV < E 1 < 
350 MeV. An enhancement from the decay tp(2S) — * rjj/tp with r\ — ► 77 in a region of 180 to 400 MeV is estimated 
by MC method, and subtracted from the data, according to the measured total number of i\)(2S) events and known 
branching fractions. The smooth background under the signal photon lines can be described by a threshold function 
from the Mri-Fit package 0] . 

bgthreshoidix) = B-(x-x ) a -e c ^-^ +c ^ x -^ 2 , (3) 

where B, xq, and a are normalization factor, threshold, and power; c\ and ci are coefficients linear and quadratic in 
x. The threshold function has been used to fit backgound with a threshold at the lower or upper side of an observed 
distribution by some experiments |19|. 

Considering the physical photon energy x = E~ t and its error y = A£" 7 due to detector resolution, the measured 
photon energy is u = x + y and its pdf function can be written as : 



h p df(u) = J c El (u - y) ■ c e f f (u - y) ■ f p df(u - y) ■ g res {y)dy = J c El (x) ■ c e //(x) • f p df(x) ■ g r es(u - x)dx, (4) 

where ce^x) — x 5 /E^ Xcj , c e //(x) = e XcJ (a;)/e(£' 7iXc/ ), f p df(x) is defined with Eq. J5J, and g r es(y) is the ICB 
resolution function. With the assumption that the El electric dipole transition for ip(2S) — > 7X c j decays (J=0,l,2) 
dominates, an E^ energy dependence is included in the folded signal shape. The detection efficiency e XcJ (a;) and 
energy resolution as a function of photon energy are included in the fitting. Normalization factors E~ UXc l and e(i? 7iXc J ) 
are photon energy corresponding to fitted Xcj mass and efficiency at the photon energy, respectively. Notice that 
parameters, masses M XcJ and widths T Xc7 (J=0,l,2), are implicitly contained in the f P df{x) function, and detector 
resolution and tail parameters are in the g res (y) function. The likelihood function, Lk(u; M XcJ ,T XcJ ), is constructed 
with three Xcj signals plus threshold background: 

2 

Lk(u; M XcJ ,T XcJ ) = bg thres hoid{u) + Aj ■ h pdf . XcJ (u; M XcJ7 T XcJ ). (5) 

j=o 

Here Aj is the observed area in each Xcj signal. 

An input-output test is performed to verify the accuracy of the fitting algorithm for the 2D problem using MC 
events. The energy dependencies of the photon detection efficiency and resolution are included in the fitting procedure. 
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A sample of MC events for the ip(2S) — > 7Xcj decays with non-zero width for both the ip(2S) and Xcj are produced. 
The photon energy distribution is fitted with the 2D pdf function convoluted with the ICB resolution. The resulting 
masses and widths of the Xcj states in this test are consistent with the MC input parameters. 

A combined fit of the three photon spectra corresponding to the ip(2S) — > 7XcO ; 7Xci> 1Xc2 decays is performed to 
three 2D pdf functions (see Eq. (J5J), each convoluted with its ICB resolution function, plus threshold background. 
The Xc2 width is fixed in the fit due to the limited statistics. The % 2 and D.F. (degrees of freedom) from the fit are 
521.0 and (520-13), where total number of free paramaters is 13. The effect of the beam energy spread |2y| in the 
measurement is also included, but is found to be negligible due to the narrow width of the i/j(2S) state. A study 
shows that the bin size (0.5 or 0.2 MeV) in the binned fits slightly affects the fitted masses and widths of the XcJ 
states. The difference in the results due to different bin sizes are added to the systematic error. The results of the 
binned fit (0.5 MeV/bin) are shown in Fig. |2 

VII. SYSTEMATIC ERRORS 

Samples of QED radiative two photon events with one photon converting to an e + e~ pair are selected for both data 
and MC simulation. The two photons are required to be emitted back-to-back. The fitted photon energy in data is 
different from the expected MC value by — 1.2ct (cr = 0.86 MeV), which has a relative error at the same level as a 
correction factor s = 0.9975 ± 0.0007 for the magnetic field [2l]]. Thus a relative precision of 0.0007 is added as the 
systematic error in the photon energy determination. 

The selection of ir° — ► 77 events with both photons converting to e + e~ pairs from 58 x 10 6 J/tjj events yields a 
data sample of 503 tt° mesons. A MC sample of 7r° mesons is generated with the same momentum and polar angular 
distributions as found from the ir° data sample. The tt° mass resolutions determined from the data and MC are in 
good agreement; their difference is 0.15±0.62 MeV/c 2 . We assume that the photon energy resolution and uncertainty 
in the direction of the photon momentum each contribute half in the 7r° mass resolution. Hence, the difference Acta/ 
of the 7T° mass resolutions between MC and data from the uncertainty of the photon energy resolution lies within 
(-0.29, +0.59) MeV/c 2 with a probability of 68.3%. We assume Aa E ^{Xc.j)/crE y {XcJ, MC) = Aa M ^ /a M o(MC), 
where ct_b t (XcJ, MC) and Acrg^ (xcj) are the MC photon energy resolution and the difference between MC and data 
for Xcj final states, and (7 M o (MC) and Actm are ir° mass resolution in MC and the difference between MC and data. 
Thus Ict confidence intervals of Act_e 7 (Xcj) for the ip(2S) — > 7XcO, 7Xci decays are estimated to be (—0.20, +0.40) and 
(—0.14, +0.27) MeV, which are used to estimate systematic errors in the determination of the XcO and Xci widths. 

The effect of the background shape uncertainty is studied using ip(^S) data and ip(2S) — > anything MC [22| . 
The relative differences in background shape parameters between floated and fixed widths of the Xc0,ci states are 
determined in fits for MC data, and fed back to correct background parameters in the fits for data. The difference 
between results for ip(2S) data with the background shape floated and fixed is taken as a systematic error. In addition, 
our MC study with non-zero width of both ip(2S) and Xcj shows that differences in the fitted masses from input values 
for the XcO and Xcl are 0.12 + 0.06 and 0.06 + 0.03 MeV/c 2 , while that for the X c2 is as large as 0.31 + 0.06 MeV/c 2 . The 
differences are attributed to uncertainties in the energy loss correction for low momentum electrons. The systematic 
errors, including the contributions from the uncertainties of the photon detection efficiency, are summarized in the 
Table D 

TABLE I: Summary of systematic errors in the determination of the \cJ masses and widths (in MeV/c 2 ). 



source 


M Xc0 




M Xcl 




M Xa2 


background shape 


0.03 


0.8 


0.02 


0.07 


0.04 


correction in magnetic field 


0.19 




0.13 




0.09 


MC simulation in as y 




+0.29 
-0.76 




+0.25 
-0.77 




different bin size 


0.02 


0.02 


0.01 


0.01 


0.02 


photon energy correction 


0.18 




0.09 




0.37 


efficiency uncertainty 


0.04 


0.03 


0.01 


0.00 


0.04 


error of M^s) 


0.034 




0.034 




0.034 


total 


0.27 


+0.85 

-1.10 


0.16 


+0.26 

-0.77 


0.39 
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VIII. RESULTS AND DISCUSSION 

With good energy resolution for low energy photons obtained using photon conversions, the precise measurement of 
the masses and widths of Xcj (J — 0, 1, 2) states from inclusive ip(2S) radiative decays can be obtained. The masses 
and widths are determined to be M Xc0 = 3414.21±0.39±0.27, M Xcl = 3510.30±0.14±0.16, M Xc2 = 3555.70±0.59±0.39 
MeV/c 2 , T Xc0 = 12.6±i:|±?:? and T Xcl = 1.39±g;fg±g:|f MeV/c 2 . The mass splittings in the Xcj(lP) triplet and their 
ratio are found to be AM 21 = M Xc2 - M Xcl = 45.40 ± 0.61 ± 0.42 MeV/c 2 , AAfi = M Xcl - M Xc0 = 96.09 ± 0.41 ± 
0.31 MeV/c 2 and p(x c ) = AM 2 i/AM w = 0.472 ± 0.006 ± 0.004. For the first time, the spin-averaged 3 Pj mass 
(weighted with the factors 2 J + 1) for the Xcj states is precisely measured in one experiment and determined to be 
M( 3 P cog ) — 3524.85 ±0.32 ±0.30 MeV/c 2 . The first errors in the results are statistical and the second are systematic. 
Correlations are taken into account in estimations of both statistical and systematic errors for the AM21, AM10, 
p(Xc) and M( 3 P cog ). Correlation coefficients between mass parameters for statistical error are obtained from the 
error matrix of the combined fit, and that for systematic error are all assumed to be 1. 

The Xcj masses (J=0,l,2) determined here are consistent with the recent measurements by CLEO JTjJ| , but have 
smaller systematic errors. Theprecisions for the Xco and Xd masses are compatible with those of previous measure- 
ments by E835 01 an d E760 while that for the Xc2 mass is not as good as theirs due to low statistics. Note that 
our Xco mass is lower than that measured by the E835 via Xco iJ/i 1 decay by 1.2 MeV (corresponding to 1.8c), but 
agrees with their later measurement via Xco — * 7r°7r° decay. The width of the Xcj states (J=0,1) determined here are 
also consistent with their values; larger errors in our widths are caused by limited statistics for both signal photons 
and inclusive 7r° mesons. 
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